Calibrating passive scalar transport in shear-flow turbulence 
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The turbulent diffusivity tensor is determined for linear shear flow turbulence using numerical simulations. 
For moderately strong shear, the diagonal components are found to increase quadratically with Peclet and 
Reynolds numbers below about 10 and then become constant. The diffusivity tensor is found to have compo- 
nents proportional to the symmetric and antisymmetric parts of the velocity gradient matrix, as well as products 
of these. All components decrease with the wave number of the mean field in a Lorentzian fashion. The compo- 
nents of the diffusivity tensor are found not to depend significantly on the presence of helicity in the turbulence. 
The signs of the leading terms in the expression for the diffusion tensor are found to be in good agreement with 
estimates based on a simple closure assumption. 

PACS numbers: PACS Numbers : 47.27.tb, 47.27.ek, 95.30.Lz 



I. INTRODUCTION 

In a turbulent flow, chemicals tend to be mixed more ef- 
fectively than in the absence of turbulence. Indeed, turbu- 
lence disperses chemicals by advecting particles along chaotic 
trajectories. This rapidly causes large concentration gradi- 
ents that speed up their mixing down toward the smallest 
scales. Turbulent mixing is a complicated and rich process; 
see Ref. (Ql] for a comprehensive review on this subject. The 
mathematical treatment of the description of turbulent mixing 
is closely related to that of turbulence itself, but it is in many 
ways much simpler and provides therefore an ideal tool for 
making conceptual progress in that field 101 • 

Here we are mainly interested in cases where it is mean- 
ingful to define a mean concentration whose scale of variation 
is large compared with the scale of the energy-carrying ed- 
dies. In such cases it can be useful to describe the change in 
the mean concentration by an effective turbulent diffusion ten- 
sor. On smaller scales the change in the mean concentration 
can still be described in such a way, but in that case the mul- 
tiplication with a turbulent diffusivity must be replaced with 
a convolution. The turbulent diffusion tensor quantifies the 
effective exchange of chemicals or other passive scalar quan- 
tities advected by the flow. If there is a gradient in the mean 
concentration C of chemicals, there will be a net mean flux 
^ = Mc of chemicals resulting from a systematic correlation 
of fluctuations in the concentration c and the turbulent velocity 
u. Here, overbars denote averaging. Under isotropic condi- 
tions with sufficient scale separation, this mean flux will be 
down the gradient of concentration, with 

^ = -AttVC, (1) 

where Kt is the turbulent diffusivity. However, modifications 
are expected when the turbulence is anisotropic. In that case 
this relation takes the form 

y^ = -««jV,C, (2) 

where is now the turbulent diffusion tensor. In this paper 
we are interested in the anisotropy caused by the presence of 
shear. One of the results one expects to see is a suppression of 



turbulent transport in the cross-stream direction. This effect is 
discussed in various physical circumstances such as geophys- 
ical flows Jit], turbulent plasmas |4], and solar physics |5, 6]. 

Much of this research is done using analytical techniques 
such as the first-order smoothing approximation and the renor- 
malization group analysis. However, in recent years it has 
become possible to calculate turbulent transport coefficients 
using numerical realizations of turbulence from direct simula- 
tions. Turbulent transport coefficients can then be determined 
by imposing a gradient in the passive scalar concentration and 
measuring the resulting concentration fluxes |7]. By imposing 
gradients in three different directions it is possible to assemble 
all components of the turbulent diffusion tensor. 

In recent years such a technique has been applied to the case 
of magnetic fields whose evolution is controlled not just by 
turbulent magnetic diffusion, but also by non-diffusive contri- 
butions known as the a effect |8, 9]. In this way it has been 
possible to investigate numerically the effects of shear and ro- 
tation in regimes that cannot be treated analytically. The tech- 
nique is known under the name test-field method, which refers 
to the fact that this approach involves the analysis of correla- 
tions for a set of different pre-determined test fields. In the 
analogous case of passive scalars, this method is now often 
referred to as test-scalar method fw\. 

Using this method, it has recently been possible to deter- 
mine the turbulent diffusion tensor in cases where the turbu- 
lence is anisotropic owing to the presence of either rotation 
or an imposed magnetic field ifToll . In the case of rotation the 
angular velocity vector ft provides a new element for con- 
structing an anisotropic rank-2 tensor of the form 1 1 li] 

Kij — kqSij + Kn^ijk^k + Knn^i^], (3) 

where Jl = r2/|r2| is the unit vector along the rotation axis 
and kq, Kf2, and Knn are functions of the flow parameters. 
Note that 17 is a pseudo vector while is a proper tensor, 
so all three coefficients in Eq. ^ are proper scalars. In the 
case of a shear flow, an obvious possible ansatz is obtained 
by replacing ft with the vorticity W = V x J7, which is 
also a pseudovector (or axial vector), and U is the mean shear 
flow. However, such an ansatz would be incomplete, because 
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it only captures the antisymmetric part of the velocity gradient 
matrix Uij, where a comma denotes partial differentiation. 
A more natural approach would therefore be to invoke both 
symmetric and antisymmetric parts of the velocity gradient 
matrix by writing it as Uij ~ Sij 



Aij , where 



(4) 



relations above to compute the five unknown coefficients in 
Eq. (|6]l via 

KsS* = K21 + Ki2, ka-S = K21 - K12, (17) 



KSsS'^/2 K22 + Kii - 2k33, 



(18) 



The latter can also be written as Ay = 
rank-2 tensor can then be expressed as 



(5) 



ei-jkWk- A proper 



= K22 - Kii, 



Kt = K33. 



(19) 



(20) 



= K^S.j + KsS^j + ni^k,, + Kss (S S)„- + KAS (A S)y , (6) Note that combinations such as nsS and kss5V2 have still 



where Kt, ks, ka, kss^ and kas are proper scalars that are 
again functions of the flow parameters. In the absence of he- 
licity, no further rank-2 tensors can be constructed from a lin- 
ear shear flow. We return to the case with helicity in Sec. lIIIDl 
An important goal of this work is to determine the coeffi- 
cients in Eq. (|6]l for a linear shear flow of the form 



U = (0,S'a;,0), 



(7) 



where S = const is the shear rate, which is not to be confused 
with the tensor S. For a linear shear flow given by Eq. O, the 
tensors S and A are constants, and their only non-vanishing 
components are 

^xy ^yx ^xy — ^yx — S/'2. (8) 

Note also that 

= -A2 = (5/2)2 diag (1,1,0), 

AS = -S A = (5/2)2 diag (-1,1,0). ^^q-. 

With these preparations we can now express all nine compo- 
nents of Kij in terms of the five coefficients in Eq. ^ as fol- 
lows: 



Kll = Kt + 4S' (ksS - KAs), 



K22 = Kt + jS (kss + KAS), 



(11) 



(12) 



the same dimension as , so in the following we shall quote 
these combinations in that form. 

In principle it is possible to construct using also the ve- 
locity vector U itself. However, U varies in x and vanishes at 
X — 0. On the other hand, we expect the components of Kij 
not to depend explicitly on position, making a construction 
in terms of U less favorable. Furthermore, the tensor UiUj, 
which has only one component in the yy position, can already 
be constructed from — A S = diag (0, 2, 0), so no new in- 
formation would be added. However, this changes when we 
also admit helical turbulent flows, because then there could 
be tensors of the form WiUj and WjUi which have compo- 
nents in the yz and zy directions. For this reason we shall also 
investigate helical turbulence in some cases. 

A comment regarding the case of rotation without shear is 
here in order. In hindsight it might have been more natural 
to write Eq. ^ in terms of the antisymmetric matrix Ay — 
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Kij — K^. Sij + Kfl^Aij + K^^(A )ij, (21) 

with coefficients that are related to those in Eq. Q via 

K^^Ko + Knn, K^ = -2Kn, kaa = 4Kon. (22) 

Evidently, this representation is equivalent to that of Eq. (O. 

In the rest of this paper we continue with the case of a pure 
shear flow. The aim is to determine the coefficients in Eq. ^ 
as functions of flow parameters such as the Peclet number and 
the shear parameter 



K33 = Kt, 



(13) 



II. SIMULATIONS 



K12 = hS{Ks - ka), 



K21 = hSins + ka), 



Kl3 = K31 = K23 = K32 = 0. 



(14) 



(15) 



(16) 



Given that all nine components of Kij can be determined from 
simulation data using the test-scalar method, we can use the 



We simulate turbulence by solving the compressible hydro- 
dynamic equations with an imposed random forcing term and 
an isothermal equation of state, so that the pressure p is related 
to p via p = pc^, where Cg is the isothermal sound speed. 
We consider a periodic Cartesian domain of size i"^. In the 
presence of shear the hydrodynamic equations for p and the 
departure U from the imposed shear flow U take the form. 



Dlnp 
Dt 



V (7, 



(23) 
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^ = -SU,y-clVlnp+f 



(24) 



where D/Di = d/dt + {U + (7) • V is the advective deriva- 
tive with respect to the full velocity, -Fvisc = p^^V • 2pi^S 
is the viscous force, i/ is the kinematic viscosity, Sy = 
^{Uij + Uj,i) — ^Sij'V ■ U is the traceless rate of strain ten- 
sor of the departure from the shear flow, and / is a random 
forcing function consisting of plane transversal waves with 
random wave vectors k such that |fc| lies in a band around 
a given forcing wave number kf. The vector k changes ran- 
domly from one timestep to the next, so f is S correlated in 
time. We have carried out simulations with helical and non- 
helical forcings using the modified forcing function 



(nohcl) 



with R,i = 



(25) 



where the non-helical forcing function. In the fully 

helical case (a = ±1) we recover the forcing function used 
m Ref. fl^l, and in the non-helical case (a = 0) this forc- 
ing function becomes equivalent to that used in Ref. ifisll . 
The forcing amplitude is chosen such that the Mach number. 
Ma ~ Unns/cg, is about 0.1. We use triply-periodic boundary 
conditions, except that the x direction is shearing-periodic. 



U{-^L, y, 2, t) = Uai, y + LSt, z, t), 



(26) 



where L is the side length of the cubic domain. This condi- 
tion is routinely used in numerical studies of shear flows in 
Cartesian geometry lfl4l[l5l] . 

In this paper we are interested in the turbulent mixing of a 
passive scalar concentration C. Its evolution is governed by 
the equation 



dt 



= -V • (UC) + kV^C, 



(27) 



where k is the microscopic (molecular) passive scalar diffu- 
sivity. In the absence of any sources, the dynamics of C 
depends essentially on initial conditions. For example, if C 
is initially concentrated in a plane with its normal pointing 
in one of the three coordinate directions, turbulence tends to 
spread this initial distribution away from the plane - regard- 
less of its orientation. Only the speed of spreading will be 
different in the different directions. The spreading is then best 
described by introducing planar averages over the same direc- 
tions as the initial distribution. These averages are denoted 
by overbars and they depend only on time and the direction 
normal to the plane of averaging, i.e. C = C{xj , t), where Xj 
denotes x, y, or z for j — 1, 3, just depending on the ini- 
tial distribution. This allows us then to quantity the speed of 
spreading by the different components of the diffusion tensor 
Kij in Eq. We do this by introducing different 'test scalars' 
and calculating the evolution for each case separately.. 

In the following we are interested in the fluxes of the pas- 
sive scalar concentration, ^ = uc, where c = C — C is the 
fluctuation around the mean concentration and u = U ~ U 



is the velocity fluctuation around the mean flow U . The test- 
scalar equation is obtained by subtracting the averaged pas- 
sive scalar equation from the original one and applying it to a 
predetermined set of six different mean fields. 



C = Co cos kxi, C =Cosinfcxi, 



(28) 



where Cq is a normalization factor. Again, the overbars de- 
note planar averaging over the directions that are perpendic- 
ular to the direction in which the mean field varies. For each 
test field C we obtain a separate evolution equation for the 
corresponding fluctuating component c'"'. 



ddPt 



= -V-([/cP9+mC^VmcP«-m^) + kV2cP«, (29) 



where p = I, 3, and q — c ov s. In this way, we calcu- 
late six different fluxes, JF^"^ = itc?"?, and compute the nine 
relevant components of Kij, 



(cos kxjF', 



sin kxjJ-'i'^) /k, 



(30) 



for i, j = 1, 3. Here, angular brackets denote volume aver- 
ages. A visualization of c^*, c^*, and c^" on the periphery of 
the computational domain is shown in Fig. [T] after about one 
turnover time for a run with k/kf = 0.1, which is smaller than 
in most of the runs analyzed in this paper. This ratio is cho- 
sen here for visualization purposes only, because this way the 
large-scale modulation compared with the scale of the turbu- 
lence becomes evident. 

We emphasize that Eq. ( |29] ) is an inhomogeneous equation 
in d"i. The term uC can be regarded as a forcing term 
that guarantees that the direction of the turbulent concentra- 
tion flux will not change with time. 

In this paper we present the values of in non- 
dimensional form by normalizing with 



s/3fcf 



(31) 



which is the expected value for large values of Pe. Here we 
have defined the root-mean-square value of the velocity fluc- 
tuation as Uinis — {u^)^/"^. 

Our simulations are characterized by two important non- 
dimensional control parameters, the shear parameter Sh and 
the Peclet number Pe, defined as 



S/{u^ 



Pe = u,.^s/{Kk{). 



(32) 



In addition, there is the Schmidt number Sc — vj k, but we 
keep it equal to unity in all cases reported below. Note also 
that in most cases we use negative values of S, so we have 
Sh < 0. The smallest wave number that fits into the compu- 
tational domain is ki — 2-11/ L. In most of the cases reported 
below we choose the forcing wave number to be 3 times larger, 
i.e. kf/ki — 3. 

The simulations have been carried out using the PENCIL 
Code |34] which is a high-order finite-difference code (sixth 
order in space and third order in time) for solving the com- 
pressible hydrodynamic equations. The test-scalar equations 
where already implemented into the public-domain code, but 
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FIG. 1: (Color online) Visualization of c^", c?" , and c''^ on the periphery of the computational domain after about one turnover time for a 
run with fe/fcf = 0.1. In the middle panel, arrows indicate the direction of the shear flow with negative S, i.e. AUy/Ax < 0. Note the clear 
sinusoidal modulation in the x, y, and z directions for the three panels, respectively. In the middle panel this modulation is already smeared 
out by the shear. 



have now been generalized to determining all nine compo- 
nents of Kij . The numerical resolution used in the simulations 
depends on the Peclet number and reaches 128"^ meshpoints 
for runs with Pe « 120. In this paper we restrict ourselves to 
time spans short enough so that the so-called vorticity dynamo 
has no time do develop; see Refs. [il6i, il7i] for details on this 
effect. 



III. RESULTS 

A. Dependence on the shear parameter 

We begin by discussing the dependence of the coefficients 
in Eq. Q on the shear parameter Sh. The result is shown in 
Fig- El for Pe — 25. It turns out that all five coefficients are 
positive. We find that Kt / ^to = const — 2 for non-helical 
turbulence and 3 for helical turbulence, independent of the 
value of shear, provided |Sh| < 0.5. The other coefficients 
show the following approximate scaling behavior: 

Ks^/Kto~5|Sh|, Kss^72Kto~30Sh2, (33) 

K;,S/Kto « 10 |Sh|3, KAs^V2Kto ~ 40 |Sh|3. (34) 

The fact that ka and kas scale with the third power of Sh 
suggests that these are higher order effects that are not easily 
captured by perturbative approaches. 

A comment regarding the values of Sh is here in order. Al- 
though values of Sh larger than unity have not yet been ex- 
plored, it is unlikely that the uprise of Kt continues. Further- 
more, one might speculate that all coefficients in Eqs. (l33T l 
and ( l34l i should eventually decrease as |Sh| — > oo. 

In Fig. |2] we have also shown results for cases where the 
forcing function has maximum helicity. No significant depen- 
dence can be seen, except for Kt which is slightly enhanced 
in the helical case with weak shear. This suggests that this 
dependence is not connected with the presence of shear. 



B. Dependence on Peclet number 

We have performed simulations for different values of the 
Peclet number and have determined the coefficients in Eq. ^ 
for each simulation. The results are shown in Fig.[3]for fixed 
Sh = 0.2. It turns out that the first four coefficients can well 
be approximated by simple algebraic functions, 

Kt ^ 2KshPe^ KsS ^ ^shPe^ 
Kto ~ Peg + Pe2 ' Kto ~ (Pe2 + Pe2)3/2 ' 

Kss-g^ _ KshPe*^ ka5 _ KshPe"^ 
2Kto "(Peg + Pe2)2' Kto " (Pe2+Pe2)2.4' 

where Ksh = 0.95Kto and Peo = 3.8 are fit parameters. In 
the case of kas the error bars are so large that no conclusive 
statements can be made. Likewise, the error bar on the first 
data point is quite large too. This is caused by the numerical 
time step becoming rather short at large diffusivities, so the 
run is short and the statistics poor. 

In Fig. m we show the dependence of the diagonal com- 
ponents of Kij on Pe. Over the range of parameters shown 
here, the difference between the three components is small, 
although there is a tendency for Kyy to be somewhat enhanced 
around Pe = 20, while Kzz is slightly smaller than Kxx- 

C. Wavenumber dependence 

We consider now the dependence of the diagonal com- 
ponents of Kij on the wave number k of the test scalar in 
Eq. ( l28T l. A dependence of Kij on k reflects the fact that there 
is poor scale separation, i.e. k/k[ is no longer small. In such a 
case, the multiplication with a turbulent diffusivity in Eqs. ([T]| 
and Q must be replaced by a convolution with an integral 
kernel |10]. In Fourier space the convolution coiTesponds to 
a multiplication. The full integral kernel can be assembled 
by determining the full k dependence and then Fourier trans- 
forming back into real space. 
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FIG. 2: Dependence of the coefficients in Eq. ((SJ on Sh for Pe = 25. 
The dashed hne in the first panel is for a run with maximum helicity. 
All runs with helicity are marked with open symbols. Filled symbols 
indicate runs without helicity. Solid lines represent the fits given by 
Eqs. l|33) and l|34j. 

The resulting dependence on fc is shown in Fig. |5] for two 
values of the shear parameter and Pe around 50. In agreement 
with earlier findings, the components of Kij show a Lorentzian 
dependence on k, i.e. 

''^^ ~ l + (aA:/fcf)2' ^^^^ 

where a « 0.2 for the k\\ and K22 components, and a w 0.4 
for the K33 component. Here, tif^ is the value for /c = 0, 
which is approximately equal to Kto^ defined in Eq. OTI ). 

Given that the Schmidt number is always kept equal to 
unity, there will be a fully developed cascade in the passive 
scalar concentration when the Peclet number is large. The va- 
lidity of Eq. (|37| i has only been tested for values of Pe up to 
60. It is unclear whether this equation holds also for large 



FIG. 3: Dependence of the coefficients in Eq. ^ on Pe for Sh = 
—0.2. The symbols give the numerical results and the solid lines 
represent fits given by Eqs. l l35t and l l36i l. 
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FIG. 4: Dependence of the diagonal components of Kij on Pe. 



values of Pe when contributions from the high wave number 
dynamics may become important in the mixing of the mean 
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FIG. 5: Dependence of the diagonal components of mhj on k for 
Sh = -0.13 at Pe = 40 (upper panel) and Sh = -0.20 at Pe = 60 
(lower panel). 

concentration. 

The case of high wave numbers is interesting in view of 
possible applications of our results to subgrid scale modeling 
in large-eddy simulations of turbulence. The highest possi- 
ble wave number is the Nyquist wave number, kj^y — tt/6x, 
where Sx is the mesh scale. In the Smagorinsky model |fl8|] 
the subgrid scale viscosity is proportional to the modulus of 
the rate of strain tensor times Sx'^. For a turbulent flow where 
the local velocity difference 6u( over a distance £ is propor- 
tional to we expect the subgrid scale viscosity to be ef- 
fectively proportional to i'^^^, suggesting an asymptotic fc"'*/^ 
scaling for k ^ fcf. Here we have identified £ with Sx and 
thus k with fcNy ■ Only for a smooth velocity field, where Su^ 
scales linearly with the separation £, the subgrid scale viscos- 
ity would be proportional to justifying an asymptotic 
scaling. This uncertainty warrants further studies of the valid- 
ity ofEq. (l32llforfc > kf. 

D. Effects of helicity 

As discussed in the Introduction, the presence of helicity al- 
lows one in principle to construct proper tensors proportional 
to WiUj and WjUi, because we have now access to a pseu- 
doscalar given by the kinetic helicity of the turbulence. If this 
does indeed have an effect, one would expect finite yz and zy 
components. In Fig.|6]we present results for Ky^ and K^y us- 
ing Pe = 25. We see that Kyz = Hzy = within error bars, so 



IV. EXPECTATIONS FROM THE r APPROXIMATION 

Passive scalar transport is closely related to the transport 
of a mean magnetic field. Commonly applied techniques 
for computing turbulent transport coefficients in mean-field 
electrodynamics are the first order smoothing approximation 
ifTgl I20I1 and the r approximation 1I 21I - I23I1 . The t approxi- 
mation consists in writing down an evolution equation for the 
quadratic correlations which, in the case of mean-field elec- 
trodynamics, is the mean electromotive force S. Its solution 
gives then an expression for S in terms of the mean magnetic 
field and its derivatives. For a recent review see Ref. \2M- 
This technique has also been used to compute the Reynolds 
and Maxwell stress in rotating shear flows ll25l - l27ll . In the 
present case of passive scalar transport one starts with the 
evolution equation for the mean flux T ~ uc, as is done in 
Refs. |i7|,i28J. Thus, we write 

—-^u.c + UiC, (38) 
ot 

where dots denote time derivatives that are given essentially 
by Eqs. (l24l i and ( |29] |. This results in quadratic and triple cor- 
relations. The sum of all triple correlations is substituted by 
a damping term of the form —Tjr on the right-hand side of 
the evolution equation for T. Here, r = St/unnsfcf is the 
turnover time and St is a positive dimensionless parameter of 
order unity (referred to as Strouhal number). This is a closure 
assumption that cannot be motivated rigorously I29tl . but it 
has been found numerically that the triple-correlations are in- 
deed locally and temporally proportional to the negative flux 
term divided by r; see Ref. |7] for passive scalar diffusion and 
Ref. [3 Oil for the case of mean-field electrodynamics. 
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As a first orientation, and in order to gain some understand- 
ing of our numerical results, we make the additional assump- 
tion that we can subsume the effects of the pressure term in 
our closure assumption. Since our forcing function f is S cor- 
related in time we have fc = and thus obtain 



— ^SSi2SikUkC + triple correlations, 



UiC = —UiUj V jC + triple correlations. 



(39) 



(40) 



The triple correlation terms result from the nonlinearities in 
the evolution equations, Eqs. (|24] | and (|29] l. In the r ap- 
proximation one substitutes the sum of the triple correlations 
by quadratic correlations, i.e. in the present case by —Uic/r 
II21II3I1I . We write the resulting equation in matrix form. 



= -UfeJ^fe - TUiUj VjG, 



(41) 



where L.^fc = Stk 
and obtain 



-(L 



ST6i26ik- We solve this equation for 



ijUk VfeC - 



dt 



(42) 



where (L^^)ifc = 5ik — Sh(5i2<^i/c with Sh = St. In the pres- 
ence of shear, the Reynolds stress tensor ujok is no longer 
diagonal, but it has finite xy and yx components. Also the 
three diagonal components are no longer the same. In the fol- 
lowing we represent UjU^ in the form 



1 








(43) 



where 5 — —UxUy/u^ characterizes the value of the off- 
diagonal components, while e — Uy/u^ ~ 1 and = 

u1/u^ — 1 characterize the change in the two lower diago- 
nal components. The dependence of 5 and e on Sh is shown in 
Fig.|7] while is found to be small. Inserting this expression 
into Eq. (l42l i. we obtain 



1 



-5 





— = ( -5 - Sh 1 + e + (5Sh 



(44) 



In the stationary state we may ignore the time derivative and 
recover Eq. (O with 



= -Sh, 



KSS^72 



5Sh. 



(45) 



(46) 



We recall that Sh is negative, and that 5 changes sign with Sh. 
Therefore we expect ks and ka to be positive, which agrees 
with the simulations. Furthermore, we expect Kyy to be en- 
hanced, which also agrees with the simulations. However, the 
slight suppression of n^z cannot be explained by the simple 
theory, because is small and perhaps even positive, sug- 
gesting at best an opposite trend. 
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FIG. 7: Dependence of 5 and e on Sh for nonhelical turbulence (solid 
symbols) and helical turbulence (open symbols). The solid line in the 
second panel has a slope of 5. 



V. CONCLUSIONS 

The present work has shown that shear introduces 
anisotropics in the diffusivity tensor for passive scalar dif- 
fusion. These additional components are proportional to the 
even and odd parts of the velocity gradient tensor, as well as 
products of these tensors. Those components that are con- 
nected with the antisymmetric part of the velocity gradient 
tensor scale with the third power of the shear parameter, sug- 
gesting that these effects cannot be captured perturbatively. 

Given that Sc = 1 in all our runs, we always have Re = Pe, 
which is at most about 100, so the inertial range of the turbu- 
lence is not very big yet. It is therefore important to investigate 
the dependence of the various transport coefficients on the val- 
ues of Re and Pe, as was done in Fig. [3] The results available 
so far suggest that the first three coefficients {nt, ks, and kss) 
do not change with Re for Re > 10. If there were indica- 
tions that the resulting coefficients change beyond Re = 100, 
it would be important to make an effort to increase the values 
of Re even further This would require more resolution and is 
obviously expensive. In view of the constancy of the first three 
coefficients, this may not be well justified. The fourth coeffi- 
cient (ka) seems to tend to zero, and the fifth one (kas) shows 
large error bars. The situation regarding these last two coeffi- 
cients may not improve significantly towards larger Reynolds 
numbers, unless the simulations are run for long enough time. 

In general, turbulent transport tends to be enhanced in the 
direction of the shear, i.e. nyy tends to be larger than k,xx and 
Kzz- Furthermore, Kzz tends to be suppressed relative to k^x- 
This is a result that is not reproduced by a simple analytical 
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closure in which triple correlations are being replaced with 
quadratic ones. In particular, there is no evidence for a sup- 
pression of turbulent transport in the cross-stream or x direc- 
tion. Instead, there is a suppression in the spanwise direction 
out of the plane of the shear flow. 

We recall that the moduli of the diagonal components of the 
turbulent diffusivity tensor are found to decrease with increas- 
ing wave number of the mean concentration in a Lorentzian 
fashion. This is in agreement with earlier findings both in 
the contexts of mean-field electrodynamics with and without 
shear ll32l [33ll . as well as passive scalar transport in the ab- 
sence of shear [10]. The limit of high wave numbers may 
be of interest for subgrid scale modeling in large-eddy sim- 
ulations of turbulence. However, it still needs to be clarified 
whether the effective diffusivity is proportional to the inverse 
Nyquist wave number to the second power, as suggested by 
our current results, or to some smaller power, ^ fc^'*/'^, as 
expected for Kolmogorov turbulence. In order to address this 
question, simulations at larger Peclet and Reynolds numbers 



are required. Such simulations do not require the presence of 
shear This is however beyond the scope of the present paper. 

Finally, we note that, in shear flows, the passive scalar 
transport properties are not affected by the presence of he- 
licity. In other words, there is no evidence for the existence 
of components to the turbulent diffusivity tensor Hij that are 
proportional to WiUj and WjUi. 
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